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We develop a kinetic theory of microcavity polaritons in presence of both Coulomb and polariton- 
phonon interaction, obeying particle number conservation. We study the growth of a macroscopic 
population of condensed particles in the lowest polariton state, under steady-state incoherent ex- 
citation of higher energy states. The collective excitation spectrum, resulting from the Coulomb 
Hamiltonian treated within the Hartree-Fock-Bogolubov framework, strongly influences the polari- 
ton condensation kinetics. In particular, for values of the excitation intensity above the condensation 
threshold, scattering from the condensate into the collective excitation modes results in strong quan- 
tum fluctuations that deplete the condensate. A numerical evaluation based on a few-level scheme 
shows that the condensate fraction is expected to be lower than 1 even far above threshold. With 
increasing system size, the role of the polariton quantum fluctuations becomes dominant, eventually 
preventing condensation to occur for system size larger than 100 /an. 

PACS numbers: 71.36.+c,71.35.Lk,42.65.-k 



I. INTRODUCTION 

Quantum fluids are the most remarkable manifesta- 
tion of quantum mechanics at the macroscopic scale. Su- 
perconductivity, superfluiditjii and more recently Bosc- 
Einstcin condensation (BEC) of diluted atoms^ are all 
examples of a system in which many particles share the 
same quantum mechanical wave function. A long sought 
and never observed quantum fluid is the BEC of exci- 
tons in semiconductors^ Presently, it is not well under- 
stood why excitonic BEC eludes experimental observa- 
tion. Presumably however, three factors are believed to 
play against BEC. First the structural disorder, which 
induces a fragmentation of the condensate, effectively in- 
creasing the condensation critical density as pointed out 
by Nozieres^ Second, the high rate of exciton-exciton 
Coulomb scattering, expected to cause a strong dephas- 
ing of the condensate already at moderate density;*^ 
which is not predicted by the standard mean-field ap- 
proach to BECiSiS Third, the strong composite boson 
nature of excitons for which, contrarily to atoms, the 
Mott transition density is rather close to the typical den- 
sities for which BEC is expected to occur iAS On the other 
hand, the possibility of achieving a quantum fluid in a 
solid-state device, with ease of control and integration, 
would open a new promising way to the implementation 
of quantum information technology^ 

Recently, it was suggested that a quantum phase tran- 
sition of polaritons in a semiconductor microcavity un- 
der steady-state incoherent optical pumping might oc- 
cur, with formation of a collective state of many polari- 
tons. M^JJ^Ll^MJMML The interest of this system 
resides in the mixed nature of polaritons, which are a lin- 
ear superposition of photon and exciton states Due 
to the energy-momentum dispersion of a planar micro- 
cavity, which is parabolic around zero momentum, the 
resulting polariton quasiparticle has a very light effec- 
tive mass at the band bottom, five orders of magnitude 



smaller than the free electron mass. Another key feature 
of microcavity polaritons is the very short radiative life- 
time, spanning the range of 1 to 10 ps depending on the 
quality factor of the microcavity. Because of this short 
lifetime and of the much slower energy-relaxation mech- 
anisms^ the polariton system under high energy optical 
excitation is strongly out of thermodynamic equilibrium. 
In this situation, the simple picture of an equilibrium 
Bose gas is completely inadequate. Incidentally, this also 
implies that the often advocated role of the light polari- 
ton mass in determining a high BEC critical tempera- 
ture, based on the simple equilibrium expression for T c , 
is actually irrelevant in all experimental situations. The 
light mass does however play a very important role in 
three other respects. First, it produces a very long po- 
lariton coherence length which averages out the struc- 
tural disorder of the semiconductor heterostructure^ 
thus eliminating the effect of condensate fragmentation. 
Second, it reduces the final-state phase space available 
for polariton-polariton scattering processes, resulting in 
a strong suppression of the dephasing compared to the 
exciton systemjSSiS! which should play in favour of BEC. 
Third, for an analogous reason it is at the origin of a 
very slow condensation kinetics^ which, in the strongly 
noncquilibrium condition of polaritons, could play con- 
siderably against BEC. This complex scenario suggests 
that polariton BEC might be possible provided that the 
density, which is required for the condensation kinetics 
to be faster than the radiative lifetime, is low enough. 
It also shows that a reliable theoretical model for polari- 
ton BEC must include polariton-polariton scattering and 
nonequilibrium kinetics on equal grounds. 

A parallel between a microcavity polaritons collective 
quantum state and conventional BEC is made hazardous 
by the Hohenberg-Mermin- Wagner theorem^S^S stating 
that a phase transition with an evident symmetry break- 
ing is forbidden in a 2-dimcnsional system. For this rea- 
son, the phenomenon has been rather interpreted as a po- 
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lariton laser transition— J-i In more practical terms, the 
BEC scenario is recovered in two dimensions if a finite 
system size is considered. In this case, condensation oc- 
curs because of the finite energy gap, due to energy quan- 
tization, between the lowest and the first excited state of 
the systemi^i^ This gap quenches the long-wavelength 
quantum fluctuations which in the limit of infinite size 
destroy the condensate. For polaritons, a finite quantiza- 
tion size is naturally introduced by polariton localization 
over a few tens of /im resulting from defects in the micro- 
cavity structure^ or by the finite size of the laser excita- 
tion spot. A few experimental results suggest some kind 
of stimulated phenomenor*i2ii&i2ii£iiLM or even a phase 
transition with spontaneous phase buildup^ but many 
observed features, among which the unexpected obser- 
vation of a thermal-type two-photon correlation function 
far above threshold^ still do not match neither the laser 
nor the BEC picture— 

Laussy et al— have pointed out that an important role 
is played by the particle number conservation. Indeed, 
in any symmetry breaking approach, a state with a well 
defined quantum phase cannot be stationary, due to the 
fluctuations of the particle number— Therefore, as in the 
theory of BEC in diluted atoms, a number-conserving ap- 
proach is needed in order to correctly describe the quan- 
tum phase diffusion of the condensate— i2LS To inves- 
tigate the appearance of condensation (either at finite 
temperature or in a non-equilibrium regime), it is im- 
portant to remark that in the BEC modelsiiSii both the 
condensed and the non-condensed phases, having differ- 
ent fluctuation terms, are considered. In these models, 
the interactions are the key feature at the origin of the 
collective excitation spectrum, and are responsible for 
the scattering kinetics that determines the ratio between 
condensate and excitation populationsiSS^S On the the- 
oretical side, many existing works on polaritons prefer 
overlooking this aspect, pursuing a strict analogy with 
the laser theoryji2*i*i or however neglecting the role of 
many-body interactions— 

In this work we develop a model of the polariton dy- 
namics which includes the polariton-polariton Coulomb 
interaction and the polariton-phonon scattering on equal 
grounds, considering a non-equilibrium steady-state opti- 
cal pump populating the high energy states. The model is 
derived within the Hartree-Fock-Bogoliubov (HFB) ap- 
proximation, as in the case of quantum fluids at finite 
temperature^*^ but imposing the particle number con- 
servation— We introduce the key assumption that 
the energy-relaxation processes are much slower than the 
polariton-field dynamics induced by the Coulomb inter- 
action. This adiabatic assumption allows to compute the 
relaxation kinetics onto a quasistationary HFB spectrum. 
This situation, which for polaritons is justified by the 
very slow relaxation kinetics within the steep region of 
the lower polariton dispersion curve, is exactly opposite 
to the case of BEC of diluted alkali atoms; 28 i 44 i 45 i 46 where 
the weak mutual interactions ensure a slow field dynam- 
ics compared to the thermalization processes. We show 



that the Coulomb interaction is responsible for a deple- 
tion of the polariton condensate even far above threshold, 
in favour of the excitations. This condensate depletion is 
strictly related to the existence of collective Bogolubov 
modes and becomes the dominant process in the limit of 
large system size, effectively preventing condensation. In 
Section II we derive the full theoretical formalism and dis- 
cuss its implications. Section III presents an application 
to a simplified few-level model, which might describe a 
situation with polariton lateral confinement and sizeable 
energy quantization. Section IV is devoted to the discus- 
sion of the numerical results. In Section V we present 
our conclusions and outlook. 



II. THEORY 

We consider the polariton in the lower branch of the 
dispersion as a quasi-particlc in two dimensions, de- 
scribed by the Bose field pk- 

\pk,P k >]= S k k'- (1) 

The lower polariton Hamiltonian in presence of Coulomb 
and polariton-phonon scattering ia 4 ^ 4 ^ 

H = h ^kp\p k + hLJ i b l b i + H c + H ph , (2) 

k q 

H c = l^2 v kk'Pl+ q pl- q Pk'Pk, (3) 

kk' q 

Hph = T,9kt(bl + b-MPk'+plPk), (4) 

kk' q 

where the the matrix element for polariton-polariton in- 
teraction Vf??, can be derived from the Coulomb interac- 
tion between excitons and from the oscillator strength 
saturation term originating from Pauli exclusion— ^ 
The polariton-phonon matrix element g^}, can be derived 
from the deformation potential interaction with acoustic 
phonons^ 4 ^ which is expected to dominate at low tem- 
perature, but could also include other electron-phonon 
coupling mechanisms. Since we aim at a kinetic descrip- 
tion of the polariton dynamics, we adopt the number- 
conserving HFB approximation. In fact, although the 
total number of particles is expected to vary in presence 
of a pump and of finite escape probability through the 
mirrors, we still cannot lift the constraint of particle- 
number conservation, basically for two reasons. First, the 
description of the condensate as a classical field would re- 
sult in an unphysical kinetic equation for the condensate, 
in which the spontaneous in-scattering term vanishes, as 
shown at the end of this section. Second, at any fixed 
time the number of particles is well defined in the real 
system. The energy eigenvalues of the Bogolubov-like ex- 
cited states depend self-consistently on the actual num- 
ber of condensed and non-condensed particles* 4 . 3 ' 46 When 
including phonon relaxation, this dependence also affects 
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the energy relaxation rates. One way to overcome the 
first of the two problems within a symmetry-breaking ap- 
proach consists in writing a separate semi-classical Boltz- 
mann equation for the condensate density and introduc- 
ing a stimulation term scaling as the inverse of the area, 
as was done by Doan et al£& Here, however, we prefer to 
adopt the number-conserving formalism which directly 
leads to fully self-consistent kinetic equations. 

In the number-conserving HFB, the polariton field is 
expressed as 



pk = PkCL+pk, 



(5) 



i.e. the sum of a condensate part P k a and a single- 
particle excitation part . The condensed particle 
operator obeys the Bose commutation rule [fi, fit] = 1 
and N c = (fit ci) defines the population of condensed 
particles, while represents the normalized wave func- 
tion of the condensate in momentum space. In a non- 
uniform condensate^, this wave function is determined 
self-consistently by imposing the relation (<Ppk) = 0, re- 
sulting in the Gross-Pitacvskii equation in the thermody- 
namic limit. The single-particle excitation p k is orthog- 
onal to the wave function of the condensate, 



kPk 



0, 



and obeys the modified Bose commutation relation 



\Pk,p{. 



Jkk' 



(6) 



(7) 



required to preserve the Bose commutation relation for 
the total field. Using these definitions, the total popula- 
tion of particles with momentum k is 



N k = (P k Pk 



\Pk\ 2 N r 



(8) 



where Nk = (p\pk) is the non-condensed population. 

The time evolution of the populations can be evaluated 
by means of the Heisenberg equations of motion. As 
a first step, we consider only the Coulomb interaction 
Hamiltonian, neglecting the polariton-phonon scattering. 
We then obtain the following equations for the dynamics 
of the field operators: 



la 



and 



k 



E 



Pk(uJkPk + vJ2ph-qPk'Pk- q ) (9) 
k' ,q 



ip k = uJkPk + v ^ Pk'- q Pk'Pk-q ~ iPka - iPfcfi. (10) 

k' ,q 



kk' 



hv, and 

we adopt mean-field factorizations, as detailed in the Ap- 
pendix. Let us introduce the amplitude of the scattering 
process bringing two particles from the non-condensate 
to the condensate 



m k .k' 



(fitfitp fc p fc , 



(11) 



Then the Heisenberg equations result in the following 
time evolution of the condensate population 



2ulm 



E Pk( & % 

k.k' .q 



Pk'Pk-q) 



2vlm< ^2 



PkPk' 



mk',k-q 



.(12) 



fc.fc' 



In order to simplify the present analysis, we specialize 
the model to a spatially homogeneous system. In such 
a limit, the condensate wave function can be safely as- 
sumed as a homogeneous function in the spatial domain. 
We expect the total deviation from this approximation 
to be small, as the actual wave-function will differ only 
at the system boundaries. In momentum space, this as- 
sumption implies P k = e 4 ^(5o,fc, where cf> is the conden- 
sate macroscopic phase. The assumption therefore im- 
plies that the condensate state is always characterized by 
k = 0. In the following the macroscopic phase factor e 4 *, 
together with its time dependence, will be included in the 
definition of the operator fi. We point out that the as- 
sumption of a homogeneous system is not in contrast with 
that of finite size, provided the system size is large enough 
to neglect boundary effects. The state orthogonality in 
this case implies p k = pk for k ^ 0. As a result, only 
the diagonal scattering amplitudes rh k = (a^ a)pkP-k) 
appear in the equations, in analogy with the anomalous 
correlations entering the standard HFB approach^ In 
this limit, Eq. ^ becomes 



Nr. 



2vlm ] rrik^j 



2ulm {to} . 



(13) 



where we have defined the total scattering amplitude m 
as a sum over all possible final states of the Coulomb 
scattering process. 

Turning to the kinetics of the non-condensate degrees 
of freedom, we define the destruction operator of a con- 
densate excitation^ 



A 1 -t~ 

A fe = -7=0 'p fc 



N 



(14) 



that creates a condensate particle by destroying a non- 
condensate one. This operator obeys quasi-Bose commu- 
tation rules 



Afc,A, 







Notice that here and in the following we assume a con- and 
tact polariton-polariton interaction, i.e. v 



A fe ,At 



= 5 k , q (N c -N k -l)/N. 



(15) 



(16) 



Introducing the standard Bogolubov transformation, the 
single-particle excitations can be expressed as 



A A 



Uk& k + V\al 



fc"-fei 



(17) 
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where U k and V* k are modal functions, and a k are 
the operators for Bose normal modes corresponding to 
the collective excitations of the systemiSi2i In partic- 
ular, the commutation rule 115|) impose the condition 
U k V k * = U- k Vl k M^2. Using Eq. 0, we obtain a direct 

relation between the one-particle density matrix (p k Pk') 
and the populations of Bogolubov modes N k = (a\d k ): 

(PlPk') ~5 kk ,{(\U k \ 2 + \V k \ 2 )N k + \V k \ 2 }. (18) 

This brings to the result, expected by symmetry ar- 
guments, that for a spatially homogeneous system the 
off-diagonal density matrix terms of the non-condensed 
states are vanishing within the mean-field approach, i.e. 



Nk,k' = {plpk 1 ) = N k 5 kik >. 



(19) 



Hence, we obtain the following equation for the popula- 
tion of the excited states 

N k = -2vlm{m k - ^(PlP-gPkP-k)}- (20) 



Finally the scattering amplitude m k obeys the equation 
(see the Appendix for the derivation) 

fh k = m k m k -2ivN k m-iv(l + 2N k )N c (N c -l) 
+ iv(l + 2N c )Y l tf$-£kf-k). (21) 

where fl k = -2[u k + v(N c - N k - 5/2)]. 

The residual two-particle correlations for excited par- 
ticles, appearing in Eqs. ij2TJ|) and (|2*TJl . require special 
care as they are the source of the off-diagonal long-range 
correlations characterizing a Bosc-Einstein condensate in 
real space. Once again, these terms in the number- 
conserving approach are analogous to the anomalous cor- 
relations appearing in a standard HFB formalism^ In 
particular, they cannot be factored in the single-particle 
basis without affecting the spatial correlation properties 
of the condensate. We will see in the next section how 
these terms can be treated within a simplified few-level 
model. Notice that the quantity rh denotes a scatter- 
ing process which would not conserve energy in a single- 
particle picture, as it describes the scattering of two par- 
ticles from one state to another one at larger energy. This 
process is actually present in our formalism because the 
spectrum of the system is modified by the interactions 
and the new eigenstates are the collective Bogolubov ex- 
citations, describing condensate fluctuations with large 
wavelength, rather than the single-particle states^ 

We now introduce the contribution to the population 
kinetics due to polariton-phonon scattering. In the limit 
of a spatially homogeneous system, we can write the 
phonon contribution to the population kinetics of the 
condensate as 



N, 



ph 



--Y 



g$((bl + b- q )(tfp k -pla)). (22) 



The Heisenbcrg equations for the operators produce 
again a hierarchy of equations for phonon-assisted corre- 
lations of all orders, which are coupled to the HFB vari- 
ables. The equations for the first order phonon-assisted 
correlations, like e.g. (b q a^p k ), are formally solved within 
the self-consistent Markov approximation. In particular, 
the higher order phonon assisted correlations entering the 
kinetic equations are factored according to the mean field 
approximation. For example 



(a) a) w 



b\p k )~(N c -l)(b\tfpk) 



(23) 



(b- q ,b\pl,p k ) ~ 5 q . q/ 8 k ,k'N k {\ + n q ), (24) 

where n q is the phonon distribution at wave vector q. In 
this way, we obtain equations of the following kind 

idt{b\o)pk) = \w k + v{N c + Nk-l)-w q \{b\a)p k ) 



(?) 

Hk.O 



(n q -Nk)N c -(l + n q )N k ,(25) 



whose formal solution can be plugged into the HFB equa- 
tions. The resulting contribution to the dynamics of the 
condensate population is: 



N, 



ph 



2TrJ2H"k+v(N c + N k -l)-w q ) (26) 

q,k 

9 k q f 



[(N k -n q )N c + (l+n q )N k 



,(fe) 



)N c + (l + n^)N k 



where gi and are respectively the effective phonon 
scattering matrix clement and the phonon population at 
the wave vector defined by momentum and energy con- 
servation. In particular, the in-plane wave vector com- 
ponent g|| is fixed by momentum conservation, while the 
z-component q z is selected by energy conservation. In 
this way, the polariton-phonon coupling introduces ef- 
fective phonon-mediated polariton-polariton interaction 
terms' the lowest-order ones being proportional to 
Ifffcfc' I 2 - T ne phonon populations n q are assumed to be 
thermally distributed at the lattice temperature. Eq. 
(|26|l is the standard Boltzmann equation expected for the 
energy-relaxation kinetics. A similar Boltzmann equa- 
tion holds for the populations N k of the excited single- 
particle states. The phonon-assisted correlations entering 
the equations for the scattering amplitudes rh kl on the 
other hand, cannot be solved analytically as was done for 
Eq. (|25|l . because Coulomb interaction couples different 
values of the momentum. 

We conclude this section by giving the explicit expres- 
sion of the phonon coupling term, as it would appear in 
the equation for the condensate if a symmetry-breaking 
approach was adopted. Expressing the polariton field as 



Pk = Pk +Pk, 



(27) 
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where P k = (p k ) is a classical field, while p k describes 
fluctuations, and considering again the uniform limit 
Pk = Sk.oP, the phonon contribution to the condensate 
field equation reads 



P 



ph 



=4 *-«)»>• 



(28) 



(j.k 



Following the same procedure adopted above in order to 
evaluate the phonon assisted terms within the Markov 
approximation, we obtain for the condensate population 
N c = \P\ 2 the following expression 



N 



ph 



= 27T2_^S(uJ k - W q ) 

q,k 



(N k -n q )N c , (29) 



Where the in-scattering term present in Eq. (|26J) is here 
instead missing. This term, responsible for the sponta- 
neous scattering into the condensate, is of course absent 
within a description of the condensate as a classical field. 
As suggested at the beginning of this section, the adop- 
tion of the number-conserving formalism prevents the oc- 
currence of this unphysical behavior. 



III. FEW-LEVEL MODEL 

As seen in the previous Section, the solution of the 
whole set of equations for the populations and correla- 
tions is a challenging task, basically because of the off- 
diagonal coupling in the phonon assisted correlations. 
In addition, the two-point correlations between single- 
particle excitations appearing in Eqs. I|2(J|I and l|21|) still 
need to be addressed in a consistent way. In a typical 
photolumincscence experiment under nonresonant exci- 
tation, the steep polariton dispersion results in a relax- 
ation bottleneck^Si with polariton population piling up 
at the boundary of the flat exciton-like region of the po- 
lariton dispersion. From there, polaritons relax to the 
band-bottom, where the actual phase transition can take 
place before they recombine emitting a photon. Given 
the very slow relaxation rates^ and the very fast radia- 
tive recombination rates, the relaxation from the bottle- 
neck region to the polariton band bottom is very likely 
to take place in one, or at most a few relaxation steps. 
Multiple relaxation steps within the steep region of the 
polariton dispersion are however very unlikely. Their 
contribution might quantitatively affect the total out- 
scattering rate from the bottleneck region, but it will 
not affect in a sizeable way the in-scattering rate in the 
ground level, which is the relevant process for final-state 
stimulation. This allows the introduction of an effective- 
level scheme, which should describe the relevant dynam- 
ics of the problem. In this scheme, the flat bottleneck re- 
gion of the dispersion, the excited states within the steep 
region and the ground state are described as three effec- 
tive levels, accounting for the respective density of states 



in the real system. We will refer to these levels as bot- 
tleneck, single-particle excitations and condensate respec- 
tively. For the reasons illustrated above, we can assume 
that the phonon mediated relaxation occurs only between 
the bottleneck level and the two other levels, while the 
Coulomb amplitudes rhk are non- vanishing only for states 
energetically close to the condensate, i.e. in the band 
bottom region. The first assumption, as already stated, 
stems from the fast radiative rate of the polariton levels 
in the strong coupling region of the dispersion, ensuring 
a negligible contribution to the in-scattering rate in the 
ground level^* The second assumption is justified by the 
fact that the Coulomb scattering amplitudes rhk can be 
important only for small wave vector k, because they are 
linked to condensate fluctuation of large wavelength, as 
discussed previously. This scheme is analogous to the one 
commonly adopted for the description of BEC kinetics of 
a diluted alkali gas 4 s . The idea behind it is that in the 
lower energy region of the spectrum, the Bogolubov field 
dynamics is the dominant process, while the low density 
of states makes the relaxation kinetics negligible. The 
opposite occurs in the higher-energy region, where most 
of the relaxation kinetics takes place but the collective 
Bogolubov excitations coincide with the single-particle 
states^ Our simplified few-level model is sketched in 
Fig. n The creation operator for the bottleneck state is 
defined by p\ (and N± is the population per mode at bot- 
tleneck), while p' now indicates the creator for a single- 
particle excitation. The corresponding total population 
of single-particle excitations is given by TV = J2kNk- 
Our purpose is to write a set of equations giving the time 
evolution of the condensate population N c , of the single- 
particle excitations population TV and of the bottleneck 
population TVi. 

As in previous treatments^ the introduction of ef- 
fective levels requires the introduction of renormalized 
coupling constants. In our case, the effective phonon cou- 
pling rates involving the bottleneck level are proportional 
to the total number of states p x in this region of the dis- 
persion. This latter is estimated from the assumption of 
a thermalized polariton distribution at the bottleneck^ 
and is related to the exciton mass M exc and the exciton 
energy thermal broadening E ksT by 



Px = {A/2n){M exc E/n 2 ), 



(30) 



(i) 



resulting in an effective phonon coupling rate g k = Px9l 
where g k is the phonon coupling rate for a scattering 
process between a state at the bottleneck and a polariton 
state with momentum k, defined as in equation i12(il) 
We now address the problem of the two-particle corre- 
lation terms entering equations (|2U|I and (|21|) . Using 
again Bogolubov transformations l|17(l , linking the single- 
particle field to collective excitations, we can write 



q . fe 



T(N) 



J2u k V k *(l + 2N k 
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FIG. 1: Schematic lower polariton dispersion. The relaxation 
from the bottleneck region to the ground state is described by 
a few-level model in which a bottleneck level, a single-particle 
excited level and the ground level are introduced. 



$> Xfe iV fe (xkN k + 2\V k \ 

k 



2\V k 



(31) 



where T(N) = N 2 [(N C + 1)(N C + 2)}- 1 and X k = 
& + 2\V k \ 2 . The quantity & = (N c - N k - l)/N 
comes from the normalization of the Bogolubov factors 
\U k \ 2 — \V k \ 2 = £,k~ One of the key points of the present 
approach consists in evaluating the Bogolubov coeffi- 
cients self-consistently as a function of the instantaneous 
populations at each time in the kinetic. This amounts to 
replacing in Eq. I|31|). the expression 



where 



it/ |2 t [Ek - (Uk + gfc)] 
of. - [E k - (uj k + <Tfe) 



E k =[(u k + a k ) 2 -4] 1/2 



(32) 



(33) 



is the Bogolubov quasi-particle energy and a k — vN£ k . 
The physical interpretation of this approach is that the 
relaxation kinetics is much slower than the Bogolubov 
field dynamics. We are therefore assuming a quasi- 
stationary situation where, at each time during the evo- 
lution of the system, a Bogolubov spectrum can be eval- 
uated from the instantaneous distribution of the quasi- 
particle populations. This is of course only allowed in 
the limit of validity of the number-conserving Bogolubov 
approach, namely for N c > N k & 

The terms in Ea. (|31|l still contain k-dependent factors 
but we can see that, for N > 1, the dominant contribu- 
tion is proportional to N 2 , via the relation 



Nk= (I + N c )N k - \V k \ 2 Xk 1 - 



(34) 



The amplitude of this term depends on the anomalous co- 
efficients Vfe , which tend to zero for vanishing condensate 
density. Therefore, in the thermodynamic limit, since the 
quantities N c and V k vanishes as expected for a Bose sys- 
tem in two dimensions, the quantity in Eq. l|31|l reduces to 
a term proportional to nN, where n is the particle den- 
sity. On the other hand, in order to solve the effective- 
level equations presented below, we want rewrite this con- 
tribution in a more useful k-independent form. To this 



purpose, we propose the following argument based on the 
two-point spatial correlation function. In real space, we 
can write 

^2(plpl q p k p- k ) = J rfrds(p t (r)p t (r)p(s)p(s)). (35) 

q,k 

Now, for a non-condensed system, at the lowest order in 
the total density, the four-point spatial correlation can 
be safely factored in terms of the two-point correlation 
(pt(r)p(s)} as 

(pt(r)pt(r)p(s)£(s)) a 2(p t (r)p(s)) 2 -(5(r-s)iV(r). (36) 
Therefore, defining the parameter 

a = N~ 2 J drds(pi(r)p(s)) 2 , (37) 



we obtain 



q p k p_ k ) ~ (2aN - 1)N, 



(38) 



In an equilibrium situation, for a non-condensed phase, 
the two-point spatial correlation vanishes at distances 



larger than the thermal length Xt = y 2irh 2 / (inksT), 
which derives directly from the equilibrium Bosc-Einstein 
distribution of excitations. In a non-equilibrium situa- 
tion like the present one, the distribution of excitations 
generally differs from the equilibrium one, but again it 
corresponds to a characteristic length. We however ex- 
pect the correlation length n of the one-body density 
function, defining the spatial correlation length for a non- 
condensed system, to be in the same range as Xt ■ There- 
fore, from Eq. H37f) we see that 



Q 



K 

~A 



(39) 



Assuming k independent of the system size, a then scales 
as A~ x £ For a system that has undergone condensa- 
tion, on the other hand, the simple factorization used 
in Ea. l|36|l is not valid, as it totally neglects the anoma- 
lous correlations. These contributions are the dominant 
term in Ea. H31|) . proportional to the products U k V k * . In- 
spection of Eq. (|31|l . shows that its dominant term still 
depends on the squared population of excitations N 2 , as 
previously discussed. We therefore propose to use the 
expression Q38H also in the condensate regime, by making 
a different assumption on the parameter a. In particu- 
lar, for a condensed system, the parameter a should no 
longer scale as in Eq. (|39|l . because of the presence of 
long-range correlations. At the same time, comparing 
equations and we see that only values a < 1 

are admitted, the limiting value a = 1 corresponding to 
a four-point spatial correlation function extending over 
the whole system size. We still expect, however, a resid- 
ual dependence of a on the system area A even in the 
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condensed regime, because by increasing the system size 
we suppress the relaxation mechanism in favour of the 
condensate depletion. Therefore, a self-consistent evalu- 
ation of the parameter a is needed for the present kinetic 
model to have the correct thermodynamic limit. We dis- 
cuss this approach below, after having introduced the 
effective-level description of the system. Concluding this 
analysis, we remark that the quantity (|38[) is real, so the 
last term in equation i|20|) gives a vanishing contribution 
when summed over k. 

In order to obtain a closed set of equations for N c , N, 
Ni and m, the energies of the single-particle excitations 
are replaced by an effective value, hui ~ Tiujk represent- 



ing a typical energy of the non-condensed levels relative 
to the ground state. This quantity plays a crucial role 
in this effective-level scheme, as it determines the finite 
energy gap that makes condensation possible in a two di- 
mensional system. We discuss below how the numerical 
value of this parameter is chosen. Similarly, we introduce 
the energy huj\ for the bottleneck states. We assume in- 
trinsic linewidths h^i for the bottleneck level and h"fo for 
the other levels, accounting for radiative recombination 
as well as nonradiative homogeneous energy broadening. 
In this way, using the relation (|38() . the resulting equa- 
tions are 



N c = -2 l0 N c + 2g c r{l + n c )N 1 (l + N c )-2g c rn c (l + N 1 )N c + 2vlm{m}, 

N = -2f N + 2gr(l + n)iVi(r? + N) - 2gTh(l + N{)N - 2vlm{m} , 
Ni = -271 JVi + 2T(g c n c N c + gnN)(l + N x ) 



(40) 



-2r 



Ni+F : 



g c (l + n c )(l + N c ) + g(l + n)(r) + N) 
-2 {2 7o + r m [g'(n' - JVi) + g"{n" - N x )\ + zft'} m 
+iv \(1 + 2N c )N(2aN - 1) - (r) + 2N)N C (N C - 1) 



where T = 71/(71 + 70), T m = 71/(71 + 370), ft' = 
[Q + v(N c + (1 — l/r))N — 5/2)] and ij is the number of 
single-particle excitation states, depending on the quan- 
tization area via the relation k = (n x x + n^y^/v^A^*^ 
We have denoted by n c , h, n' and n" the phonon popula- 
tions defined by energy conservation within the Markov 
assumption, as shown in equation l|26|) . Correspond- 
ingly, g c , g, g' and g" denote the effective phonon cou- 
pling strength at the same energies, renormalized by the 
number of bottleneck levels p x . The equation for the 
scattering amplitude contains an oscillating term, whose 
frequency depends on the actual energy needed to cre- 
ate condensate particle. Notice that all this quantities 
depend on the actual condensate and non-condensate 
densities and vary self-consistently during time evolu- 
tion. In particular, the parameter a can be obtained 
sclf-consistcntly at each step of the kinetics, by equating 
expressions (1311) and l|38() . The effective- level represen- 
tation of the condensate excitation naturally implies an 
effective Bogolubov factor V, via Eq. and corre- 

sponding values for the other related quantities. The 
expression for a takes then the compact form 

a(w) = i±f (i +2 jqfV (41) 

2 + N c \y rj x 2 J 

In particular, we find again that the first term in paren- 
thesis scales as A' 1 , consistently with the previous dis- 
cussion. Actually, it turns out that the precise value 



of a is not critical in determining the condensation dy- 
namics, provided that the two limits a = and a = 1 
are not reached. In particular, in the next section, we 
will show how the self-consistent result does not differ 
much by one one obtained using a constant value for a. 
We have introduced a steady-state pump rate F in the 
equation for the effective bottleneck level. This quantity 
represents the number of particles per unit time and per 
state, which enters the system following the relaxation 
from higher-energy states, as in a typical experiment with 
non-resonant continuous-wave excitation>i2iiii^ii^ For a 
given total pump flux / (total rate of polaritons per unit 
area), the quantity F can be rewritten as F = fAp^ 1 . 
This shows that for a given pump flux the quantity F 
does not depend on the area A of the system, as can be 
argued from Eq. lj3U|) . 

The phonon-mediated interaction results in a Boltz- 
mann dynamics^4 with scattering coefficients g c T, gT, 
g'T rn and g"T m . In the kinetics described by our model, 
the initial growth of the condensate is triggered by final- 
state stimulation within the relaxation process. It takes 
place in the ground state which is separated by an en- 
ergy gap ui from the excited level. Above the stimula- 
tion threshold, a crucial role is played by the Coulomb 
scattering amplitude m, which enters the first two equa- 
tions in 140f) with equal factors up to a sign, and de- 
termines the actual condensate growth or depletion. As 
we will see, this term determines the macroscopic con- 
densate fraction reached by the system. Indeed, we can 
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compute the approximate analytical steady state frac- 
tion, neglecting the oscillating term in the equation for 
rh and making the high density limit N 3> 1,N C 3> 1 
with Ni/N c — * 0. This is justified by the fact that, 
above threshold, Ni is expected to reach a fixed value 

(thr) —1 

Ni ~ FthJi ■ Neglecting the phonon population and 
taking T m g' ~ r m g" = g m and Tg c = g c , we obtain for 
m the steady state solution 

_ i = i v(2aN-N c )N c N_ 
2 7o - g-mFthli 
This result, plugged into equation for 7V C , gives 

N c „2aN~2 (7 ° = g£^L^9 Z 9mF thl ^) 

V 2 N 2 V / 

predicting the coexistence of condensate and non- 
condensate even far above threshold. This behaviour, 
as will be also seen in the numerical solution, originates 
exclusively from the Coulomb scattering, whereas a pure 
Boltzmann dynamics would always predict a condensate 
fraction approaching 1 above threshold. 

The adoption of this simplified model calls for some 
additional remarks, in order to understand its limita- 
tions and the relevance of the parameters. As already 
argued, we can expect that the relaxation from the bot- 
tleneck region into the low-energy states is qualitative 
well described within a few-level approximation, as the 
inclusion of other intermediate states would only result 
in a finite increase of the pump threshold. This in turn 
would not introduce significant changes to the exciton 
density at bottleneck and to the condensate density, as 
they mainly depend on the phonon coupling rates for 
the direct scattering processes. The Coulomb interaction 
results both in a modified spectrum (with subsequent 
variation of the relaxation rates) and in the occurrence 
of coherent scattering processes between the condensate 
and the states close to zero momentumiS Both effects 
are mainly important for the lowest-lying states. In or- 
der to correctly describe the dominant role played by the 
low-energy states, the parameter Two has to be set to a 
value lower than a simple average of the excited energy 
eigen-values. This parameter represents in our model the 
energy gap between condensate and excited states. It 
thus plays the same role as the finite-size energy quanti- 
zation in a fully two-dimensional system of finite size. We 
follow this prescription and therefore set the parameter 
hcu to the energy gap that would result from the system 
area A. The presence of more excited levels would im- 
ply an increase of the total coherent outscattering rate 
with respect to our simplified assumption which therefore 
sets a lower bound to the effectiveness of the condensate 
depiction mechanism. 

IV. NUMERICAL SOLUTION 

For the numerical evaluation, parameters of a typi- 
cal AlGaAs microcavity with one embedded GaAs quan- 



tum well have been used, with a quantization area A — 
100 /im 2 . The quantity A enters the definition of the 
phonon coupling terms 2 ^ as well as the expression for 
and implies a number of excited levels r\ = 30jS for 
which we take the representative energy Tuo = 0.1 meV 
(see the discussion in the previous section). For these pa- 
rameters, the effective phonon coupling strengths linking 
the bottleneck and the low-lying states are g c ,g ~ 1 peV, 
resulting in a very long relaxation time (r re i ax ~ 100 ps), 
consistently with the adiabatic assumption that we have 
made. The Coulomb matrix element is evaluated to be 
hv = 5 x 10~ 4 meV— The other parameters are T = 
10 K, fi l0 = 0.2 meV,?i7i = 1 meV, Kl R = 3.5 meV, 
n(u)!-w) = 1.9mcVand p x = 10 4 -S, The value h^ = 0.2 
meV is typical for the radiative lincwidth of a polari- 
ton at zero exciton-cavity detuning in such a system. 
The value h'ji = 1 meV for the exciton-like part of 
the polariton branch, on the other hand, accounts for 
the density-dependent exciton dephasing rate^ expected 
in the vicinity of the transition density resulting from 
our calculations. We have solved numerically the set 
of equations as a function of time. For each value of 
the pump density /, we observe a time-dependent tran- 
sient followed by a stationary solution for all the quan- 
tities. Fig. [2t a ) displays the stationary populations as 
a function of /. In order to check that the choice of 
parameter a is not critical, we compare the results ob- 
tained using a constant value a = 0.4 with the results 
obtained calculating sclf-consistently a from Eq. I|41|l . 
at each step of the relaxation. The figure clearly shows 
that the two methods give very similar results. Obvi- 
ously, by using a density-independent a, the condensate 
fraction at low pump density is overestimated). Below 
threshold, Boltzmann relaxation results in an increase of 
the polariton population, in which the condensate pop- 
ulation remains microscopic, i.e. N c <C N. A threshold 
occurs at / = fth, for which N c ~ 1. Above threshold, 
the bottleneck population Ni reaches a saturation value 
while the condensate population N c becomes a macro- 
scopic fraction of the non condensed polariton popula- 
tion N, which in turn continues to grow with /. The 
behaviour of the ratio N c /N is plotted in Fig. Hfb). At 
/ fth the condensate fraction approaches a finite value 
lower than 1, consistent with a = 0.4. Fig. [2{c) displays 
the imaginary part of the steady-state Coulomb scatter- 
ing amplitude rh. Close to threshold, in correspondence 
to low values of the condensate fraction, this quantity 
takes positive values, thus favoring condensation. Above 
threshold, on the other hand, it takes large negative val- 
ues, resulting in condensate depletion. Hence, Coulomb 
interaction plays a crucial role during phase transition, 
as expected according to both laser— and BEG— *2*^ quan- 
tum theories. In order to clarify the effect of the Coulomb 
interaction, we compare in Fig. Ufa) and (b) the steady 
state solutions obtained neglecting all Coulomb terms in 
1)400. By inspection of Eqs. (|40|) . it is clear that in this 
case the populations obey a standard Boltzmann dynam- 
ics. Without Coulomb interaction, therefore, a standard 
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FIG. 2: Steady state solutions vs normalized pump intensity 
for a — 0.4 (thin lines) and for a calculated self-consistently 
(thick lines): detailed discussion in the text, (a) Solid line: 
N c (plotted also on a log- log scale in the inset). Dashed line: 
N. Dot-dashed line: Ni. The same quantities iV c (dot-dot- 
dashed) and TV (dotted), computed neglecting the Coulomb 
interaction, are also plotted, (b) The corresponding conden- 
sate fraction (solid line) compared to the result computed 
neglecting Coulomb interaction (dashed line), (c) The imag- 
inary part of the Coulomb scattering amplitude Im{m}. 



three-level Boltzmann equation is recovered. At thresh- 
old, the condensate population starts to grow due to final 
state stimulation and the system undergoes a complete 
transition to a fully condensed regime with a condensate 
fraction equal to 1. 

The threshold fth depends on the energy gap us. In 
the limit of a system of infinite size, this energy gap van- 
ishes and fth becomes infinite. In order to understand 
how the size of the system affects the condensate frac- 
tion, on the other hand, in Fig. [3]we compare the results 
obtained with different values of A. For this compari- 
son, we assume that Q scales as A' 1 , according to the 
energy quantization in a quantum system of finite size. 
The results plotted in Fig. 0show that the condensate 
fraction reached far above threshold depends on the sys- 
tem size. This condensate depletion is due to the co- 
herent scattering terms out of the condensate that orig- 
inate from the kinetic HFB equations. In the thcrmo- 
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FIG. 3: Condensate fraction as a function of the pump in- 
tensity for different values of A. A = 100 /im 2 (solid line), 
A = 200 fim 2 (dashed), A = 400 /im 2 (dot-dashed) and 
A — 600 /im 2 (dot-dot-dashed). The pump intensity is renor- 
malized to the threshold value for A — 100 /im 2 . In the in- 
set we show the asymptotic condensate fraction above pump 
threshold, as a function of the system area. 



dynamic limit, therefore, our model predicts a vanish- 
ing condensate fraction, as imposed by the Hohcnbcrg- 
Mermin- Wagner theorem^! If in the present model we 
neglect the Coulomb interaction, the coherent scattering 
terms are absent, resulting in a condensate fraction that 
always tends to 1 far above the threshold pump intensity 
(see Fig. [21 (b)). The present kinetic model therefore 
correctly reproduces the behaviour of a two-dimensional 
Bose system in the thermodynamic limit, where the long- 
wavelength fluctuations are expected to destroy the con- 
densate. 

To better understand the interplay between scattering 
amplitudes and condensate growth, we display the time 
evolution of N c and N in Fig. 0fa) and of the imagi- 
nary part of fh in Fig^Jb). At short times, the quantity 
have not a stationary value, due to the relaxation dy- 
namics and the condensate population is small. Corre- 
spondingly, the quantity Im{m} takes positive values. At 
longer times, when all the quantities reach a stationary 
value, with a macroscopic occupation of the condensate, 
Im{m} has turned to negative values, implying the re- 
sulting condensate depletion. 

The Coulomb interaction is therefore the key mecha- 
nism determining the condensate fraction under steady- 
state conditions. For the parameters used in Fig. [21 the 
saturation value of N\ corresponds here to a bottleneck 
polariton density tib ~ 2 x 10 1:L cm -2 , larger than the 
optical saturation density, suggesting that condensation 
cannot occur in a simple system with one quantum well. 
When using the parameters of the experiment by Deng et 
alfisil& i.e. considering a sample containing 12 quantum 
wells at T = AK with hflji = 7.5 meV, we obtain a qual- 
itatively similar behaviour (not shown) with Ni = 20 at 
threshold, resulting in a bottleneck polariton density per 
quantum well ~7x 10 9 cm~ 2 , lower than saturation 
density and in fairly good agreement with the experimen- 
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FIG. 4: Time evolution of iV c (solid line in (a)), N (dashed 
line in (a)) and Im{m} (in (b)) for a pump value / = 2f t h- 



a value significantly lower than 1, depending on the sys- 
tem size, as expected for a two-dimensional Bose system. 
The condensate fraction can even become vanishing if the 
size exceeds a few tens of /im. 

We conclude that the coexistence of condensate and 
non-condensate, caused by the Coulomb interaction, is a 
dominant aspect of the polariton condensation dynamics. 
This feature, unpredicted by models based on a standard 
Boltzmann kinetic approach, can affect strongly the co- 
herence properties of condensed polaritons and possibly 
prevent condensation in the most common experimen- 
tal conditions. This result holds great importance in the 
light of the numerous experimental claims of polariton 
Bose-Einstein condensationi2iiii^iiSiiSi24 and of the re- 
cent achievements in lateral confinement of microcavity 
polaritons over the micrometric scalei^iSi 

Acknowledgments 

We are grateful to I. Carusotto and R. Zimmcrmann 
for fruitful discussions. We acknowledge financial sup- 
port from the Swiss National Foundation through project 
N. 620-066060. 



tal estimatei 14 ' 1 ^ We argue that, in the experiment by 
Deng et ai, the unusually high values measured for the 
two-photon correlation function might be explained by 
the coexistence of condensed and non-condensed phases 
even far above threshold. 



V. CONCLUSIONS 

We have presented a kinetic theory of microcavity po- 
lariton condensation. The theory is based on a number- 
conserving HFB description of the polariton quantum 
field, treated using a density-matrix formalism. This al- 
lows inclusion of various scattering processes, inducing 
the energy relaxation and the final-state bosonic stim- 
ulation which causes the condensate to grow. In par- 
ticular, we describe the polariton-phonon scattering, but 
other scattering mechanisms^2iS could in principle be in- 
cluded. Differently from BEC kinetic models in atomic 
physics, here the HFB field dynamics is much faster than 
the relaxation kinetics. As a consequence, the steps of the 
relaxation kinetics follow adiabatically the spectrum of 
collective modes stemming from the HFB field equations. 
We show that in this non-equilibrium regime, the conden- 
sate fluctuations strongly influence the phase transition 
even far above threshold. In particular, they induce scat- 
tering from the condensate to the excitation modes that 
can result in a small condensate fraction in presence of 
slow energy relaxation rates. 

Within a few-level model, we have performed a nu- 
merical evaluation of the condensation kinetics. It turns 
out that in realistic experimental conditions the conden- 
sate fraction above threshold approaches asymptotically 



VI. APPENDIX 

The equations of motion appearing in this paper are 
derived in the mean-field limit, i.e. factoring higher order 
correlation terms in single-particle or Bogolubov quasi- 
particle population terms. In this appendix we give the 
details of this prescription. Let us consider the sys- 
tem having a defined number of condensed and non- 
condensed particles, N c and N = J2^k respectively. 
We compute the following expectation values of the two- 
particle quantities 

(tftfaa) ~ N C (N C - 1) 



and 



N-J2N, 



(44) 



(45) 



kk' qq' 



Consistently, we can introduce an approximated Hamil- 
tonian in terms of linearized operators, having the same 
expectation values as obtained in mean-field limit. The 
linearized operators are 



(46) 

a^aa ~ N c a, (47) 

AA'ViVi 1 - Nk,qp{,Pq' + N Kq ,p\,p q (48) 

p\pl q p h ~ [N q , h -8 q>k )pl q + N- q>k p\, (49) 

Pl+q'-kPq'Pq ~ N q+q ,- k , q .p q + N q+q ,- k , q p q , (50) 

- 8 q+q >- ktq S qtq >N q p q , 

PkPkPg — NkPq + N k>q p k - 5k, q N k p k . (51) 
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These relations can be used to factor the correlations 
entering the equation for the scattering amplitudes rh k 

fh k = -2i(u) k - luq - 2v)rh k + v(ata\afa + aa f )p k p-k) 

- v + rfajplP-qPkP-k) 

+ v{bJa)aa{j$_ k p- k + p k p\)) 

+ v^(a ] a\f>l +q ,_ k p- k +p k pl +q , +k )p q p q ,). (52) 

qq' 

In detail, we apply the following factorizations 

(a^a^a^apkp-k) — (N c - 2)rh k (53) 

(a)a)aa)p k p_ k ) ~ (N c - l)rh k (54) 

(a f a f plp q p k p- k ) ~ (N q - 8 q . k N k - 8 q - k N- k )fn k 

+ N q ,kfh q -k + Nq-kfhq^ (55) 

((aot + otoJpJptaPfcP-*) - (l + 2iV c )(ptp t _ 9 p fe p_ fc ) (56) 



(attfaa(pl k p. k + p k p\)) ~ AT C (JV C - 1)(1 + N- k + N k ) 

(57) 

(a?a*P q+ q>-kP-kPqPq') - 2N q+q ,_ k ^m q _ k 

+ N q+q >-k,-krhq,q' 

- 6q tk 6q>,kN k mk) (58) 

(a^a^ p k p t q+ql+k p q p q/ ) ~ 2N q+ q l+kt q,fh q - k 
+ N q+q i +k ^ k fh q ^ q i 
- 8- k ,q'8- k ,qN- k m k ). (59) 

Notice that we can also rewrite the double sums as 

y^^Ng+q'-fcrfrhg-k = ^ N qtk m q - k + ^ iV^/mfc 

qq' q,q'—k q' ,q—k 

- N k m k + ^2 N q+q ,_ kiq ,rh q - k 

q,q'^k 

(60) 

and the last term can be neglected, within the assump- 
tion of a spatially homogeneous system. In this way, Eq. 
(12 It is recovered. 
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